Convergence to equilibrium in competitive 
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Abstract. We study a generalized system of ODE's modeling a finite num- 
ber of biological populations in a competitive interaction. We adapt the 
techniques in [S] and [2] to prove the convergence to a unique stable equilib- 
rium. 

Resume. Nous etudions un systeme generalise d'equations different ielles 
modelisant un nombre fini de populations biologiques en interaction compe- 
titive. En adaptant les techniques de [8j et [2], nous prouvons la convergence 
vers un unique equihbre stable. 

Version frangaise abregee. 

Nous etudions le comportement en temps grand de modeles de dynamique 
de populations. On considere un nombre fini de sous-populations, correspon- 
dant chacune a un trait ou type different. Ces populations interagissent entre 
elles de fagon competitive. En notant ni{t) I'effectif de la sous-population 
numero i, un des modeles les plus classiques est le systeme de Lotka-Volterra 
competitif 

d 

-^rii = (n -'^bijnj)ni, i = l...N, 
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ou bij > 0. On se place ici dans le cadre plus general du systeme 



d , , 



y Ki{a) L\S^Bj{a)n^{t) \ dP{a) 



ni{t), 



1...N 



avec {Q, P) un espace mesurable. Ce systeme pent s'interpreter comme un 
modele avec ressources generalisees. 

En utilisant les techniques developpees dans [8] pour une version continue 
du premier modele, et dans [2], on pent facilement montrer 

Theoreme Supposons que L est une fonction surW, positive snr IR+, que 
K et B sont des fonctions positives appartenant d L°°{dP{a)) fl L^{dP{a)) 
et que 

(z) (Competition striate) L est strictement croissante et pour tout I < i < n, 

ri < J^Ki{a)L{oo)dP{a) oil L{oo) := lim^^+oo -^^(a:) G (0,+cxd]. 

(u) (Symetrie) II existe Ci > tq Bi{a) = CiKi{a) 

{Hi) (Non extinction) Pour touti, r^ > J^Ki{da)L{0) dP{a) 

(iv) (Non degenerescence) Pour I C {1 . . . A^}, soifR^ I'ensemble des n G 

tels que rij = pour tout i ^ L Pour tout Id {1 . . . A^} il y a au plus un 

neR^ tq 



ri- j Ki{a) L Bj{a) dP{a) 



0, \/i e I. 



Alors 3! n 



.riN) G 



\ {0}, tel que pour toute solution n{t) 



{ui, . . . , Tijsf) du modele generalise avec une donnee initiale ni(0) > Vi, on 



n{t) — )■ n, quand t — )■ +oo. 



En particulier ce resultat implique 

Proposition Supposons que > pour tout i et que la matrice bij verifie 
3C G (M;)^ tq Ci bij = bji Cj, et ^ u, Uj bij Q > Vii G \ {0}, 

Alors 3! n = (ni,...,nAr) G MJ^ \ {0} tel que pour toute solution n{t) = 
(ni, . . . , Un) du premier modele avec donnee initiale nj(0) > Vz, 

n{t) — > n, quand t — )■ +oo. 
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1 Introduction 



We study the long time behaviour of models of population dynamics. We 
consider a finite number of subpopulations whose dynamics is governed by a 
system of competitive ODEs (in the sense of Hirsch, see e.g. [6]). We denote 
hj Hiit) , i = 1 . . . N , the number of individuals of the subpopulation i. 

The most classical models are competitive Lotka-Volterra equations 

d 

■^Ui = {n -'^bijnj)ni, i = l...N, (1.1) 

j 

where bij > 0, and the models with a finite number of resources 

= {-di + '^Ikr]ki)ni, (1.2) 

k=l 

where r]ki > and the Ik are given by the Rolling II functional response 

ik — - 



2^i=i Vki rii 



This type of system appears in biology when one studies the dynamics of a 
system of interacting species (see [3 El |9]). It also appears in Trait Substi- 
tution Sequence models, where one considers a population structured by a 
continuous phenotype (see equation (11.41) . (11.61) on this matter), where only 
a small number of traits are present (see flU\ [1]). These models have been 
used to develop the theory of Adaptative Dynamics (see flOi [Tl [3]). 

Previous asymptotic studies on this type of equations concern either very 
general properties (the existence of a carrying simplex [6]), or precise results 
but only for low dimensional systems (A^ < 3 ^llj), under strong assumptions 
of the coefficients (for instance, the matrix (bij) is supposed to be diagonal 
dominant, see [7]), or only on local properties (the equilibrium population is 
locally stable, or populations rii do not vanish). 

Note that both equations (11.11) and (11.21) may be interpreted as discrete 
versions of continuous models. To each subpopulation corresponds a pheno- 
typic trait G M'^, and then posing 



N 



n{t,x) = ^ni{t)5^„ (1.3) 



1=1 
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one finds that Eq. f 1 1.1 1) for instance is equivalent to 



dtn{t, x) 



r X — 



b{x,y)n{t,dy))n{t,x) 



:i.4) 



r(Xi] and Oij = b{xi,Xj] 



with Tj = r(xj) and bij 

The long time behaviour of the continuous model fll.4p (with bounded 
initial data instead of Dirac masses) was studied in [8j. For a symmetric b 
defining a positive operator, the convergence to the unique stable equilibrium 
was proved. For the case with resources, the result is essentially contained 
in [2], which generalizes the derivation of [4|. 

The study of the discrete or continuous models corresponds to slighty 
different biological questions; in the continuous case, it is for instance con- 
nected to the issue of speciation, or how from a continuum of traits a few 
well separated ones (the "species") are selected; in the discrete case, one is 
rather concerned about survival or extinction of each subpopulations. From 
a rigorous mathematical point of view, a result in the continuous case does 
not imply anything for the discrete one. However it is easy to apply the 
techniques developed in |H] and [2] to the discrete models; that is our aim. 
First of all, we consider the very general equation 



i- J K,{a)L \^ B,{a) n,(t) j dP{a) 
with (fi, P) any measurable space, or in the continuous case 



ni{t), i = l...N, 
(1.5) 



dtn{t, x) 



r{x) 



I K{x,a)Ll / B{y,a)n{t,dy)) dP{a) 



n{t, x). 



:i.6) 



We prove the following 



Theorem 1 Assume that L is on R and non negative on R_|_, that K 

and B are non negative, in L°°{dP{a)) fl L}{dP{a)) and that 

{%) (Strict competition) L is strictly increasing and ri < Ki{a)L{oo) dP{a) 

for all 1 < i < N , where L(oo) := lim2,._s.+oo L{x) G (0, +oo]. 

(m) (Symmetry) There exists Cj > s.t. Bi{a) = CiKi{a) 

{iii) (Non extinction) For any i, r-j > j^Ki[a)L{Q) dP{a) 

{iv) (Non degeneracy) For any subset I C {l...A^}, let M.^ be the set of 
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n G s.t. rii = for all i ^ I. For all I G {1 . . . N}, there exists at most 
one ri G M'^ s.t. 

n- j Ki{a) L Bj{a) dP{a) =0, Vi G /. (1.7) 

Then there exists a unique n = (wi, . . . , un) G M.^ with n 0, s.t. for any 
solution n(t) = (ni, . . . , n^) to (11. 5p with initial data ni(0) > for any i, 

n{t) — n, as t —> +oo. 



Note that, if one had ri > f Ki{a)L{oo) dP{a) for some i, then nj(t) — )• +oo 
if nj(0) > 0. Assumption (z) hence ensures the non-explosion of the system. 

Eq. (II. 5p could be directly derived from simple biological considerations. 
It assumes that the reproduction rate of a population of type i (or with 
trait Xi) is the difference between a fixed rate depending only on the trait 
and a competitive interaction with the other populations, resulting from 
the interaction with the environment. The state of each component of this 
environment (indicated by different values of a) is given by the sum 

J2B,ia)n,it). 

j 

Each such component has some independent effect on the reproduction. To 
get the total reproduction rate one sums over those. 

Eq. (II. 5p is hence an obvious generalization, with a possibly infinite num- 
ber of resources, of the model (II. 2p . It also contains the Lotka-Volterra 
system (II. ip . In this case. Theorem [T] gives 

Proposition 1 Assume that > for all i and that the matrix bij satisfies 
3C G (M;)^ s.t. d bij = bji Cj, and ^ m uj C^ > G \ {0}, 

ij 

(1.8) 

then there exists a unique n = (rii, . . . ,njv) G with n ^ 0, s.t. for any 
solution n{t) = (ni, . . . , n^) to (II. ip with initial data ni(0) > for any i, 

n{t) — )■ n, as t ^ +oo. 
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This result shows that, in Lotka-Volterra systems which are symmetric 
in the sense of fll.Sp . the competition between a mutant trait and a resident 
population leads to a unique stationary state, regardless of the initial popu- 
lation state. This is precisely the assumption needed in [1] to apply a limit 
of large population and rare mutations to an individual-based model. In 
particular, Thm. 2.7 of pQ applies to symmetric competitive Lotka-Volterra 
systems. 

Proof of Prop. [H Define the matrix rriij = Ci bij. Note that m is symmetric 
and positive definite. Hence there exists an orthonormal basis of eigenvectors 
U\ i = 1 . . . N, and corresponding eigenvalues Aj > 0. 

Then put L = ld,n = {!,... ,iV}, P = ^Eii<^., = v^^^", 

Ki{a) = C^^ \^\^U" and note that 



N ^ N ^ 



a 
1 

a 



N / N \ 



j Ki{a)L Bj{a) dP{a). 



Therefore Eq. f ll.Sp indeed yields f ll.ip in that particular case. 

Conditions {€) and {ii) of Theorem [T] are obviously satisfied. Conditions 
[Hi) holds since > for all i and 1^(0) = 0. As for condition (zv), assume 
that for a subset / one has two vectors nj, 7 = 1,2, s.t. = for z ^ / and 



V2 e /. 



Put 5n = — and simply note that Srii bijSrij = pour i = 1 . . . N. 

This means that dn = and proves (iv) since Ylfj=i ^i^ij ^''^i ^^j = 0- 
Hence the proposition is implied by Theorem [TJ □ 



Note that the same argument works in the continuous case and Eq. (II. 4p 
is a particular case of f ll.Gp for x G O a bounded domain. The condition on 
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b is 



C{x) b{x,y) = C{y) b{y,x) 



C (x) b{x , y) n{x) n{y) dx dy > Vn 7^ 0. 
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One still puts L{^) = C,- Notice that C{x) b{x,y) defines a compact, selfad- 
joint and positive operator on L'^{0). Diagonalizing the operator, one gets 

C{x) b{x,y) = ^Xifaix) faiy), 



with Aq, > the eigenvalues, tending to +00 and /„ the corresponding nor- 
malized eigenvector. It is hence enough to take f2 = N and K{x, a) = 

In the particular case where b{x, y) = b{x — y) on the whole and C = 1, 
by Fourier transform, the condition on b means that b > 0. One then takes 



K{x, a) = (cos(a ■ x) y b{a), sm{a ■ x) y b{a)) 



2 Proof of Theorem [T] 

The proof is based on the study of the following Lyapunov functional 



Fin) 



/TV \ N 

Tj 77, j, 

•'^ \j=l J i=l 



where H is an antiderivative of L and hence strictly convex. 



(2.1) 



2.1 F is a strict Lyapunov functional 

Let n be a solution to (11. 5p . Then by a direct computation 



d 



N 



Hi 



1=1 



j^K,{a)L{^B 



j{a)nj ) dP{a) — rj 



• (2.2) 



Therefore F{n{t)) is non increasing and its derivative in time vanishes only 
on stationary solutions to (II. 5p . i.e. -F is a strict Lyapunov functional for the 
system (II. 5p . 



7 



Thanks to condition (z), 



1^ > a IU{a)L(^^^n^ dP{a) - r,^ > a > 



if rii is large enough. Therefore, there is a constant a' > s.t. VF{n) ■ 
n > a'\\n\\ if ||n|| is large enough. This implies that F{n) — > +oo when 
||n|| — )■ +00, and entails that n{t) is uniformly bounded. 

Let n G M.^ be a steady-state of fll.Sp and let / be the set of i s.t. Ui > 0. 
Then, for any i E I one needs to have 



J^Ua) L ^^B,{a)n^ dP{a) 



By condition (iv) there is at most one such solution for every /, and there are 
only a finite number of possible J, F has then a finite number of steady-states. 

Classical Lyapunov functionnals' techniques then entail that the solution 
n{t) to fll.Sp converges to a steady-state n for any initial condition n(0). 

2.2 The functional F is convex 

Compute 

N 



d^F 



driidrik 
Hence as L is increasing 



j^Bi{a) Bk{a)L' \y2B,{a)nA dP{a) 



(2.3) 



i,/c i \j = \ / 

F is therefore convex and any local minimum on M.^ is global. 

Since (11.51) has a finite number of stationary solutions, this clearly implies 
that F admits a unique global minimizer n. Otherwise, F would reach its 
minimum on the whole segment linking two distinct minimizers. 

The object of the next subsection is to prove that n satisfies a stronger 
property: this is the unique ESS of the system. 
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2.3 Uniqueness of the ESS 

Any local minimizer n G M.^ of the functional F necessarily satisfies 

J Ki{a) L Bj{a)nj\ dP{a) = r^, Vi s.t. rii > 0, 

^ (2.5) 



/ Ki{a) L Bj{a)nj 1 dP{a) > ri, Wi s.t. rii = 0. 

This condition corresponds to the usual definition of an Evolutionarily Stable 
Strategy in adaptive dynamics (see for instance [3]). It turns out that there 
exists at most one ESS, n. Hence being an ESS is a necessary and sufficient 
condition to be the global minimizer of F. 

Indeed take two n' G M^, 7 = 1,2 satisfying fl2.5p and compute 

> ^ a n] (^i - j Ki{a) L B,{a)n]^ dP{a) 
+ ^ Q ^2 ||r, - ^ K,{a) L B,{a)n]^ dP{ 



This last quantity is equal to (thanks to fl2.5p ) 

+ ^ C. (n^ - n]) (^r, - ^ ^a) L S,(a)n] j dP{ 

L |^fjS,(a)njj -L |^fjS,(a)n]j j dP{a). 

As L is strictly increasing, this implies that for P a.e. a, ^^i-Bj(a) {n} — 
nf) = and by (iv), it means that n-^ = n^. 



and to 
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2.4 Conclusion of the proof of Thm. [T] 



Assume that ?ij(0) > for all 1 < i < N. We know from Subsection 12.11 
that n(t) converges to a steady-state n when t oo. 

If h does not satisfy fl2.5p . there exists i E {I, . . . , N} such that 



Since nj(0) > 0, > at all times, and the linearized equation around n 
shows that n cannot converge to n: 



provided that ||n — n|| is small enough. 

Therefore, h = n, and the proof of Thm. [T]is completed. 
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